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Abstract 

Using simple and general arguments we propose an effective Hamiltonian for the 
description of low-energy pure QCD. The Hamiltonian is a function of spatially con- 
stant collective modes. Its eigenstates can be organized into bands classified by the 
irreducible representations of an 0(8) group. The latter also determine parity and 
charge conjugation of the states. The energy spectrum agrees well with the glueball 
spectrum as measured on the lattice, and in particular the level ordering with respect 
to spin is naturally explained. 
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1 Introduction 



QCD is widely considered to be the fundamental theory governing the strong interactions. 
Its chief success to date has been in describing experimental data at high energies — deep 
inelastic scattering — where the partonic nature of hadronic constituents plays a dominant 
role. In this regime the asymptotically free nature of the QCD running coupling constant, g, 
comes to the fore allowing for a perturbative expansion using as basis states free, i.e. g — > 0, 
gluons and quarks. At low energies, however, comparable succcess has been conspicuous by 
its absence. As is well known this is due to the phenomenon of infrared slavery which has as 
a consequence the confinement of coloured states, i.e. gluons and quarks are bound together 
to form hadrons. How to describe the crossover as a function of scale between such radically 
different effective degrees of freedom remains as one of the most challenging problems in 
particle physics. 

Without reliable access to the infrared via the fundamental theory itself recourse is often 
made to treatments that use additional assumptions, such as adding confinement by hand (as 
in Refs. 0-fl), or that use effective field theory descriptions such as working directly with 
the meson and baryon fields and adjusting the interaction parameters to experimental data 
@. The fundamental question remains however: how does one derive the effective theory 
consistently from the underlying more fundamental theory? In the present QCD has 

only one parameter, Aqcd, how do the interaction parameters of the effective theory depend 
on it? Finding the answer to this question is, of course, a very ambitious enterprise which 
will not be resolved in this paper. 

One notable attempt, related to the approach of this paper, is the work of Luscher, van 
Baal and collaborators || who use a "non-perturbative" approximation of Born-Oppenheimer 
type wherein one tries to integrate out "fast" (short wavelength) fluctuations thereby leav- 
ing an effective theory for the "slow" (long wavelength) modes. In the case of QCD if this 
separation is carried out considering the theory to be confined in a small box with periodic 
boundary conditions, which acts as an infrared cutoff, then the fast modes can be integrated 
out using perturbation theory. In fact, as the QCD vacuum is translationally invariant, we 
can consider integrating out all non-constant modes thereby yielding an effective field theory 
of constant gauge fields. In effect, one is now considering a one-dimensional field theory, 
i.e. quantum mechanics. In the case of the pure gauge theory the spectrum of the resulting 
Hamiltonian will lead to the glueball spectrum. However, this is the glueball spectrum associ- 
ated with a "universe" that is smaller than a fermi. To describe the real world it is necessary 
to take a box that has at least a few times the size of the glueball. Unfortunately, this leads 
to a breakdown in the perturbation theory that is used to integrate out the non-constant 
modes. The physical reason is simple: the perturbation theory is associated with free glu- 
ons, however, the long wavelength modes are better described by other non-perturbative 
effective degrees of freedom. Thus, although the methodology in principle is a very powerful 
one, leading to an effective field theory that depends only on the running coupling constant, 
in practice one finds that in trying to extend it to larger boxes one encounters the same 
fundamental problem that led to the introduction of the box in the first place. 

Here we will take a somewhat less ambitious path motivating an ansatz for an effective 
Hamiltonian via an analysis of the group-theoretical structure of QCD, then analyzing some 
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of the consequences of this ansatz, in particular the spectrum of glueball masses. First of all 
we will take the effective Hamiltonian to depend on spatially constant modes only. From the 
point of view of the functional integral in principle it is always possible to integrate out the 
non-constant modes. However, if one is implementing an approximation technique to do it 
one has to check that the resulting approximation is trustworthy. As mentioned this has been 
investigated by Liischer, van Baal and collaborators ||. Their approximation breaks down 
when long wavelength gauge modes (other than the zero mode) become strongly coupled. 
Nevertheless, given the translation invariance of the QCD vacuum one would still expect an 
expansion around the ground state to be an expansion around a spatially constant state. 
We will take that attitude here and work with an effective Hamiltonian that only depends 
on spatially constant modes. We will restrict attention to the pure gauge sector of QCD. 
Hence, when considering excited states we will only be considering glueballs. Given the 
present experimental status, it is difficult to compare our results to the outcome of actual 
measurements, however, we can compare with results obtained from lattice simulations that 
restrict to the pure gauge sector. 

As a first step, in section 2, we formulate a classification scheme for a many-gluon system, 
assuming the dominance of only one mode with J p = 1~. This will enable us to order many- 
gluon states into bands defined by the irreducible representations of an 0(8) group. This 
step alone plays an important role, allowing us to explain the level ordering observed in 
|| [T(|. This result is to a large extent model independent, the only assumption being that 
gluon number (here we are referring to constituent gluons not perturbative g gluons) is 
approximately a good quantum number, and that the more gluons in a state the higher the 
corresponding energy. 

Next, in section 3, we make an ansatz for the effective Hamiltonian of QCD in terms of 
spatially constant gauge fields motivating the use of a kinetic energy that is quadratic in 
momenta. We further introduce a change of variables on the space of spatial gauge potentials 
such that the new variables are related to the intensity of the spatial gauge field, the spatial 
quadrupole distribution and the rotational angles in space and colour space. These variables 
have been widely used in nuclear physics [[11], |12| and also in QCD || Though we 



work with constant modes here, this restriction is not necessary in general. In fact, the new 
variables as functions of space and time qualify as candidates for effective fields that better 
describe low energy QCD. With a simple ansatz for the potential term the resulting effective 
Hamiltonian turns out to be a function of a one-dimensional oscillator, a five-dimensional 
oscillator and some Casimir operators of symmetry groups appearing in the problem of many 
gluon systems. The first oscillator describes changes in the intensity of the gluonic field while 
the second is related to the quadrupole deformation of the gluonic field with respect to an 
"intrinsic" system yet to be defined (it will be similar in structure to the "body fixed frame" 
of a nucleus). One feature of the effective Hamiltonian is that there exists a conserved 
quantity which one can naturally identify with the number of constituent gluons. 

The resulting effective Hamiltonian depends on several parameters which we determine in 
section 4 by fitting the corresponding spectrum to several glueball states with spin, parity and 
charge conjugation J PC and masses known from lattice gauge calculations. We subsequently 
make predictions for all other states and compare, where possible, with lattice simulations 
and results from other phenomenological models, such as the MIT bag model. Finally in 
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section 5 we draw our conclusions. 



2 Gluon Taxonomy 

As a first step to understanding the structure of the spectrum of a many gluon system one 
has to determine which values of J PC can appear, where J refers to spin and P, C to parity 
and charge conjugation, respectively It is implicitly understood that only colour zero states 
are allowed. To deal with this problem in an efficient way it is not sufficient to consider the 
SU(3) colour and 5*0(3) spin groups only.Q One reason for this is that gluon fields do not 
transform according to the fundamental irreducible representation (irrep) of the colour group. 
The main point, however, is that we are dealing with a system of identical bosons, hence 
it is necessary to implement Bose-Einstein statistics. This can be done most conveniently 
by considering higher groups. With respect to colour rotations, the gluonic fields transform 
like vectors in an eight-dimensional space so we can associate a rotation group 0(8) with 
colour transformations. Extending the latter rotations to unitary transformations, we arrive 
at the group U (8) ||14|| . In the same way the spatial rotations can be associated with an U (3) 
group containing the 50(3) spin group. Bose-Einstein statistics are then implemented by 
connecting at this level the spin and colour transformations formally to an {7(24) group. 



We thus obtain the following group chain jl5| : 



[N] [hh 2 h 3 
17(24) 
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where for each group the quantum numbers of the corresponding irrep are given and £ and 
K are multiplicity indices appearing in the respective reductions. The symbol [N] denotes 
the completely symmetric irrep of U (24) for N gluonic modes, which implements the Bose- 
Einstein statistics. The possible irreps of {7(8) and {7(3) are then given by the Young 
tableaux [hih^h^, consisting of three rows with hi, hi and /13 boxes, respectively, where 
hi + h 2 + h 3 = N. The fact that the {7(8) and U (3) irreps are coupled to the completely 
symmetric irrep of {7(24) constrains their Young tableaux to be equal ||14|| . 

The {7(3) group reduces via SU(3)j, where J denotes spin, to the rotation group 50(3). 
The reduction rules are given in Ref. fl6| , but will be summarized here for easier reference. 
First the reduction from the irrep [hihihz] of {7(3) to (p,q) of SU(3)j is simply given by 
•p — h\ — h% and q = h 2 — h^. The possible irreps J (with magnetic quantum number M) of 
50(3) are then determined by the possible values of the multiplicity index K, 

K = min(p, q), min(p, q) — 2, . . . , or 1 , (2) 



1 Note that we use SO(3) rather than SU(2) as we are restricting our attention to integer spins. 
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and 



J = max(p, q), maxQo, q) — 2, . . . , or 1 for K = 

J = K, K + l,...,K + max(p,q) for AT ^ . (3) 

The group that describes rotations in colour space is U(8), and it reduces via 0(8) to 
the physical colour group 577(3), where we restrict to colour singlets, i.e. to the trivial 
irrep (0,0). The intermediary orthogonal group 0(8) will play a very important role in the 
following developments. The 0(8) quantum numbers are (a^u^^O) , where only three of 
the four numbers can be non-zero as a consequence of the restriction of the U (8) irrep to 



three rows [fTTH - The largest 0(8) irrep contained in a given irrep [hih 2 h^\ of U{8) is the one 
with = hk, k = 1,2,3. The precise rules for the reduction of U(8) to 0(8) are given, 
for instance, in Ref. [|17] (there is also a computer program available for arbitrary irreps 
of different classical groups []I8"| ). The reduction of 0(8) to £77(3) will be carried out in 
a recursive manner using the reduction of U(8) to 0(8) just mentioned and the reduction 
of U(8) to 577(3). The latter is given explicitly for up to five (constituent) gluons in Ref. 
||15|| , and the reduction for an arbitary irrep of U(8) can be found in Ref. fT9| . The recursive 
procedure for the reduction of 0(8) to SU(3) is described in detail in Appendix A. The results 
are given for up to six constituent gluons in Table 1, considering only those irreps of U(8) 
and 0(8) which contain at least one colour zero irrep. Finally, note that the multiplicities 
appearing in the reduction of £7(8) to 0(8) and of 0(8) to SU(3) are denoted globally by £ 
in ®. 

The above classification is not applicable in its present form to cases where more than 
one gluonic mode is relevant. It is, however, always valid when one particular - - even 
non-constant — mode is considered to be dominating at low energies, subject only to the 
condition that it has spin 1, colour (1, 1) and, in view of our further considerations, parity 

(-!)■ 

In order to discuss the properties of the states under parity and charge conjugation we 
introduce boson creation and annihilation operators, b£ and b l a , for the gluon field A^, i.e. we 
write A l a ~ (b^ + b^) with a proportionality factor depending on certain parameters. What 
we will present in this section is independent of these parameters, which will be fixed in the 
next section. The operators b£ and b l a satisfy the usual commutation relations. 

A basis of states can be characterized unambiguously by the quantum numbers in (|I|) 
(some of them are redundant, see also below). Any such state can be obtained by applying 
the pair operators 

tf/ = EM >< (4) 



[■I] 

M 

in (01) represents the spin coupling of two tensors A and B to spin J |2T1, where in our 



to a minimum weight state, i.e. a state with = u>k, k = 1,2, 3. The notation [A x B] 



case J = 0, 2 due to the spin and the bosonic nature of the gluons. Note also that the pair 
operators q^ are 0(8) scalars because of the summation over a. The conjugate operators 
q$ annihilate the minimum weight states.^ 

2 Togcther with the generators of the U(8) group these operators form the algebra of the symplectic group 
Sp(6,R), which, however, will not be considered further here (see Ref. \12 and references therein). 
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A general basis state can then be decomposed as follows: 



Af 

The symbol p(qj^™ 1 ' 2n2 ' 2n ^ KlJl represents a coupling of (n\ + n 2 +n 3 ) gluon pair operators 
from Eq. @ to the £7(8) (or £7(3)) irrep [2ni, 2n2, 2723] and to spin J\ with multiplicity index 
K\. Due to the commutativity of the pair operators the Young diagram [2ni, 2n2, 2713] must 
have an even number of boxes in each row ||17|| . The operator V(q') is subsequently coupled 
with a minimum weight state to the total £7(8) (or £7(3)) irrep [/11/12/13] and to spin J, the 
£7(3) (or equivalently £7(8)) coupling being indicated by Cg>. This coupling is such that 

N = hi + h 2 + h s = 2(m + n 2 + n 3 ) + u-y + u 2 + w 3 . (6) 

It is natural to interpret the quantum number N as the number of constituent gluons. One 
should, however, be careful not to identify N with the number of the constituent gluons 
which arise in other phemonenological models, as discussed in section 4. At any rate, Eq. 
(P) gives the precise definition of N in our model. 

We did not explicitly indicate in Eq. (|5]) the quantum numbers N and (p, q) , which are 
uniquely determined by hi,h 2 ,hs, the latter via p = h\ — h 2 and q = h 2 — h^. Furthermore 
it is understood that the state in (^) is a singlet with respect to the colour S77(3)|] The 
multiplicity label £ is replaced by £' in the second representation in (|5]) because part of 
the multiplicity is taken care of by the different couplings of [2ni, 2n 2 , 2n 3 ] with [cj!^^] to 
[hih 2 hs\. The division of the state into a coupling of {n\ +n 2 + n 3 ) gluon pairs and the "rest" 
is related to the concept of generalized seniority, where a particle system can be divided into 
pairs of particles coupled to colour zero and spin or 2 and the rest where no pairs of these 
types appear. The term "generalized" is used because seniority is normally associated with 
a coupling to a total scalar. 

The key point of the decomposition (|5]) is that the gluon pair operators change neither 
parity nor charge conjugation, which are therefore already determined by the minimum 
weight states, i.e. the ones with [/ii/12^3] = [^1.^2^3] • Furthermore, it turns out that the 
minimum weight states carry definite parity and charge conjugation, as will be demonstrated 
in the following. It is due to these facts that the 0(8) irreps play a prominent role in the 
classification of many-gluon states. 

The parity of an iV-gluon state is readily determined by noting that each gluon field A % a 
- and therefore each operator — carries parity (—1). Therefore, the parity of a state 
with 0(8) irrep (wiu^a^O) is given by 

Obtaining the charge conjugation of the states is considerably more involved. The reason 
can be seen by considering the simple example of a one-gluon state. The properties of A l a 
(and hence of b^) under charge conjugation are deduced in Ref. |^TJ and turn out to be 



CA l a = rj a A l a (no sum over a) 



3 In the pure gauge theory, the gluons have to be coupled to colour (0, 0), which is not necessarily true in 
the presence of quarks. The notation in Eq. (||) can be readily extended to the latter case. 
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with 



_ , -1 for a = 1,3,4,6,8 
+1 for a = 2, 5, 7 



in the standard representation of SU(3), where the generators are proportional to the Gell- 
Mann matrices. Since one-gluon states belong to the U(8) and 0(8) irreps [1] and (1000), 
respectively, we have an example of an 0(8) irrep lacking a definite charge conjugation. As 
we will show in the sequel one can, contrary to this observation for one-gluon states, associate 
definite charge conjugations with all colour-zero states. 

In order to see that, we have to construct the states from so-called elementary permissible 



diagrams (epds) ||22|| . An epd is, in our case, a minimum weight state as considered before, 
coupled to colour zero and definite spin, which is elementary in the sense that it cannot be 
decomposed into other couplings. It will be denoted by 

(N, [t^waws], J) = [6 f ® • ■ ■ ® tfp:p ]J , (10) 

where N = u)\ + 0J2 + uj^ and the 0(8) irrep is (0U1U2OU3O). Note that we always couple 
to M — J, which implies that a product of two epds, say (N, [^1/12^3], J)(N', [hxh' 2 h' 3 ], J'), 
represents a state with spin J+ J'. However, such a product of epds has still to be projected 
onto a definite U(3) (or equivalently U(8)) irrep, meaning that the epds are defined here 
with respect to the maximum weight of SO (3) and not of U(3) (£7(8)). 

Any minimum weight state can be represented by products of all possible epds. Including 
the pair operators of Eq. (|J) (with M = J) in the list of epds, we can then via construct 
any state with the quantum numbers of (Jl]) from the coupling of epds.[] There exist systematic 
procedures to determine the total number of epds, and the application of these general 
methods to the case at hand (including states with colour which are important in the presence 
of quarks) is presently worked out |24]]. However, for the time being, we will adopt a more 
pedestrian approach which is completely satisfactory in the case of small gluon numbers. It 
should be noticed that in general the choice of the set of epds is not unique being somewhat 
analogous to the choice of a basis for a vector space. 

The epds with up to three constituent gluons are (summation over repeated indices is 
understood) 



(2, [2], J,) = [blxbljW 
(3,[3],J 2 ) = d abc \[blxbl]^xbl 



[J* 
J 2 



(3,[1 3 ],J 3 ) = fa bc [[blxblpXbt] [JS] , (11) 
L J J3 

the first one being the pair operator (|j), while the others are minimum weight states in the 
sense defined before. The spins denoted by Jk are limited to the values J\ = 0, 2, J 2 = 1, 3 
and J3 = 0. A complete list of epds with up to six gluons, as well as the decomposition of 
some minimum weight states with seven or eight gluons to be used in the discussion of the 
glueball spectrum, can be found in Appendix B. 

4 As a detailed example for the construction of states with the help of epds see the discussion of the 
five-dimensional harmonic oscillator in Ref . |E3| . 
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We can now use the decomposition of the states into epds to determine their charge 
conjugation, using the fact that the epds are charge conjugation eigenstates. For example, 



by construction the epd of two gluons in ( |TT]) has charge conjugation C = +1. For the epds 
with three gluons we have for the symmetric coupling C = — 1 and for the antisymmetric 
coupling C = +1, which can be readily verified taking into acount the combinations of 
(a,b,c) which appear in these couplings and using the list in @. In general, as shown in 
Appendix C, the total charge conjugation is simply given by 

C={-l) nd , (12) 

where is the number of ci-symbols appearing in the coupling. Using this fact and the list 
of epds in Appendix B we can deduce the charge conjugation of all states constructed via 
Eq. (H) from minimum weight states with up to six gluons. The results are given in Table 1. 



3 An Effective Hamiltonian for QCD 

In this section we will motivate an ansatz for an effective Hamiltonian for purely gluonic 
QCD at low energies. The Hamilton density in the temporal gauge is of the form 

n = - 1 -Y,(doA:) 2 + V(A), (13) 

where the first term is the kinetic energy density and the second is the potential which 
includes everthing else, such as (V x A a ) 2 . 

To obtain the effective action of the quantum theory one must integrate over all possible 
physically inequivalent configurations of the gauge field. This is, except in certain special 
cases, impossible to do exactly. The most popular approximation technique is perturbation 
theory, however, for a discussion of QCD at low energies this simply will not do. Another 
possible approximation, as mentioned previously, is of Born-Oppenheimer type || wherein 
one tries to integrate out "fast" (short wavelength) fluctuations thereby leaving an effective 
theory for the "slow" (long wavelength) modes. Here, we make an ansatz for the low-energy 
effective Hamiltonian, taking it as a function of spatially constant gauge fields only. Given 
the translation invariance of the QCD vacuum, this is not physically unreasonable. The basic 
assumption behind this ansatz is that the dynamics of the constant modes alone can give 
us information about the glueball spectrum, i.e. the excited states of the theory, although 
obviously the inner structure of the states cannot be described by constant modes. As was 
shown in the work of Luscher, van Baal and Koller ||, this assumption is correct for the 
system located inside a box with periodic boundary conditions and a size of up to 0.7 fm 
(approximately) . 

After having integrated over all non-constant modes our general ansatz for the quantized 
effective Hamiltonian is 

H "'-^lwir + v '" {A) - (14) 
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where the effective potential V e ff is a complicated function of A % a that will be modelled by a 
simple ansatz further below. The approximation effected in (|i~4[) consists in neglecting higher 
derivatives with respect to A, as well as terms mixing derivatives with powers of A. This 
adiabatic approximation has been justified in the case of a finite volume and to one-loop 
order |J, and here we will simply assume that for our purposes the same approximation can 
also be applied to the more general case. The constant B e is a wave function renormalization, 
which we will treat as an adjustable parameter. 

After choosing the temporal gauge and integrating out all non-constant modes there 
is still a residual gauge symmetry left in Eq. (|T3D, namely the invariance under spatially 
constant, i.e. global (and time-independent) gauge transformations. We hence impose on 
the physical states the condition 

G a V phys (A) = , (15) 

where the G a are the generators of global gauge transformations. Eq. (|l"5|) is equivalent to 
the statement that physical states have no colour, in which form the condition was imposed 
on the many-gluon states considered in the previous section. The consistency of the present 
gauge fixing procedure with the integration over all non-constant modes was shown in Ref. 

We will now change variables. The transformation is chosen such that part of the new 
variables is intimately related to the gluon pairs introduced in the previous section. The 
interpretation of these variables will be given below. Explicitly the transformation is given 



by 



4 = Ee*0S(0)ASW0. ( 16 ) 

k=l 

where the constant modes of the vector potential A l a depend on the spatial indices (i — 1, 2, 3) 
and the colour indices (a = 1,2, ...,8). D fe f is the rotation matrix in the spatial three- 
dimensional space with cartesian components and Agi fe a is the corresponding rotation matrix 
in the colour space of dimension 8. Of the latter only the last three rows appear. The angles 
4> a b for rotations in colour space can be chosen in such a way that only 18 angles appear in 
Eq. (Tl6|). Together with the three angles 9ij for spatial rotations and the three Qk we then 
have (formally) 24 degrees of freedom as is the case for the A l a on the l.h.s. of Eq. (|IJ]). Eight 
of the angles (p a b correspond to global gauge transformations, hence the generators G a in Eq. 



(|i5|) can be identified with the angular momentum operators for the corresponding rotations. 
The condition ([15]) then simply shows that physical states are independent of these eight 
angles, thereby effectively reducing the number of degrees of freedom to 2 x 8, as might have 
been anticipated. 



Transformations similar to (|1J]) have been used in other areas of physics [IT], |12[ and 
also in the work of van Baal and Koller JS|. In both cases the coordinates were taken to 
be space independent as in Eq. (]T6|). However, Eq. (]T6| ) can be readily extended to the 
more general case where the coordinates include space (and time) dependence. In any case, 
the interpretation of this coordinate transformation given below is the same irrespective of 
whether the coordinates are space dependent or not. 

We will now consider the interpretation of these new variables. We first define the 
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composite field 



a 

= Y. D kMQlD\]{B). (17) 
k 

Eq. (pj]) gives the rotation in three-dimensional space of the matrix (g^) into a system in 
which the matrix is diagonal, which will be called the "intrinsic system". The matrix (g^) 
has both monopole and quadrupole components which are obtained by simply coupling to 
the appropriate spin, i.e. 

gg ] = lj | JM) Qij , (18) 

hi 

where J = 0, 2 indicates the spin and (li, lj \ JM) is a Clebsch-Gordan coefficient coupling 
two cartesian vectors to a spherical tensor of spin J. Eq. ( |T8"D represents of course the 
analogue to the gluon pair operators qffi of the previous section (see Eq. (^)). For J = 0, 

gives the square of the intensity of the vector field. The spin 2 part describes the quadrupole 
distribution of the intensity of the vector field. Thus the Qk carry the information about 
the monopole and quadrupole distribution. In Eqs. (f!T)| ) and ( |TTD the angles % give the 
orientation of the intrinsic system with respect to the laboratory system and can be defined 
in many different ways. For instance, instead of the rotation angles 8ij in two-dimensional 
planes one could just use Euler angles. The angles 4> a b describe analogous rotations in colour 
space. 

The monopole-quadrupole tensor (%•), being quadratic in the gauge field, can only de- 
scribe excitations with an even number of gluons. In addition, by the summation over colour 
indices in (|I7|), the variables and thus the excitations they describe do not change parity nor 
charge conjugation. Of course, "non-pair" excitations can still be described via ([16]) and will 
be classified by 0(8) irreps. 

In summary, our conjecture is that the variables Qk, 4> a b and Oij represent effective degrees 
of freedom that will be useful in describing the collective behaviour of the gluon field, in the 
same way as in nuclear physics analogous variables describe the collective rotational and 
vibrational modes of the nucleus. We believe that a generalization of these variables to 
include space-time dependence could also be useful in QCD. 



As described in Ref. ||12|| , instead of the variables Qk one can use g, b and c which are 
related to the former by 



2 

Qk= 3 



2 r 27T 



lb cos c k 

3 



(20) 



where k — 1, 2, 3. The parameter b gives the absolute deformation and c its deviation from 
axial symmetry. There are restrictions however [12], |25| : Firstly, the sector from c = to 7r/3 



already represents all possible physically different situations. The other sectors (six in total) 



can be mapped bijectively to the one described The Jacobian of the transformation 
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( pOD becomes singular on the boundaries between different sectors. Secondly, the restriction 
that the left-hand side of Eq. ( p0|) be positive imposes a c-dependent upper limit upon the 
variable b. 

It will turn out to be even more convenient to replace b by the variable (3 related to it 



via p| 




1 -b 2 

(21) 



- 3b 2 + 2b 3 cos(3c) 

with the scale factor a defined in Eq. (^) below. The range of (3 is from to oo (see Ref. 
|25fl). For b or {/3/a) small we have (3/ a ~ b. 

Now we will rewrite the Hamiltonian of Eq. (|14D in terms of the new variables of Eqs. 



For the time 

being we will concentrate on the kinetic term T = H e ff — V e ff. In 
effect, for this case the transformation to the new variables has already been performed in 
Ref. |l2j where the number A, referring to the number of nucleons, has to be replaced by 9 



and in particular the orthogonal group 0(A — 1) becomes 0(8), otherwise nothing changes. 

We will in the following briefly outline the procedure and give the main definitions and 
results. In the first step the kinetic term T is rewritten in terms of the variables Qk {k = 1, 2, 3) 
and the components of the spin operator in the intrinsic system associated with the angles 
9ij. T also depends on the generators C^ + k 1 ,^+k 2 (ki,k 2 = 1,2,3) which form an 5*0(3) 



subgroup of 0(8). As shown in Refs. |12| and p6| , the matrix elements of the generators 
approach for large 

a 2 = ujx + u 2 + u 3 + 6 (22) 

the matrix elements of the generators Ck {k = 1,2,3) of the 50(3)' subgroup of an U(3)' 
group in the irrep [u;iu;2u;3] . The primes indicate that a contraction procedure |27] was used 



in terms of the expansion parameter (1/cr 2 ). The [wiu^u^] irrep of U(3)' contains the irrep 
(A = iQ\ — L02, ji = uj 2 — uj-s) of its 577(3)' subgroup, and the latter can be reduced to 50(3)' 
by using the rules given in Eqs. (Q) and (^). We will identify the components of £ with 
angular momentum operators later on. 

The volume element in terms of the new variables is 

{Qi ~ q\){Qi ~ qI){q\ ~ Ql){QiQ2Q-ifdQidQ 2 dQ^dVtedVt <)) , (23) 

where dVlg and dVL^ refer to the volume element of the respective angles. Their explicit form 
depends on the specific angles chosen (Euler or others). We have also chosen Q\ > g 2 > Q3 
corresponding to the sector in the (b, c)-plane with < c < tt/3 in order to avoid ambiguities 
in the choice of the intrinsic system. 

Next, we change the volume element via the transformation 

$ = (QiQ2Q3) 5/2 (Qi + Q 2 2 + Ql) 2 * (24) 

to the new states $ in the Schrodinger picture. The change of volume element is such that 
for small (3/a ~ b the new volume element takes the form 

/3\ 4 

— I sin(3c) dgdpdcd^le dVts (25) 
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as given in Ref. [25], neglecting corrections to the next order in (/3/cr). 

Implementing these changes, the kinetic term in the effective Hamiltonian subsequently 
acquires the form 

W d 2 R 2 \ 

where R 2 is an operator which contains linear and quadratic derivatives with respect to the 
variables (3 and c. It is also a function of the spin operators in the intrinsic system and of 
the last three generators of the 0(8) group, replaced by Ck after the contraction discussed 
above. Additionally, it depends on combinations of Casimir operators of the 0(8) group and 
of its subgroups 0(7), 0(6) and 0(5) in the canonical chain. 



In Ref. [J2 an expansion is made in powers of the parameter (1/cr ) defined in Eq. (|22|). 
In the case of Ref. |12j this expansion parameter is always very small. In the case at hand, 
where it starts from 1/6 and decreases with larger u>k [k = 1, 2, 3), it is not too small, however 
one would expect the expansion not to be too bad. Actually, we will expand not only in 
powers of (1/cr 2 ) but also in powers of ((3 /a), leading to the contributions summarized by 
the formulas (5.6) to (5.11) of Ref. [Q. However, it is necessary to keep the term given by 
Eq. (5.11) in this reference (R 2 below) because it can be, in our case, of the same order as 
cr 2 . We thus arrive at the following expression for the R 2 operator of Eq. (p6[): 

R 2 = R 2 + R 2 . , (27) 



where 

.2 



R 2 « 2(T 2 



9 fl4 9 1 9 ■ ro ^ 9 ^ L l 

$ ™ - 7^7—7^77^ sm ( 3c )7i7 + E 



/3 4 0/T d(3 (3 2 sm(3c)dc v ' dc ^ 2/3 2 sin 2 (c - 2vrfc/3) 




and 



with 



(21 



R 2 w 2C 2 (A,/i)-3^/: 2 (29) 

fc=i 



Lfc — Jfc+^fc 
C 2 (A,/i) = A 2 + A/i + /i 2 + 3(A + fi) 

A = Ldi - u 2 , [i = uj 2 - u; 3 , (30) 

the ujk being the quantum numbers of the 0(8) irrep. 

The above is a good approximation in the case that (/3/cr) is very small. However, as 
this is not entirely the case we will modify the expression such that the contributions from 
higher orders in (/? / cr) are taken into account via a redefinition of the interaction parameters 
(a procedure commonly used in the collective model of nuclei PS]). At this point our model 
begins to look even more phenomenological, similar to a Landau- Ginzburg ansatz, where the 
theory contains parameters which have to be adjusted to some kind of "experiment' 
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case to lattice calculations. Nevertheless, the kinematical structure is maintained with the 
hope of learning something about the spectral structure of pure gluonic QCD. 

In order to illustrate the parameter redefinition we consider the term in square brackets 
in Eq. It can immediately be identified with the kinetic energy associated with a 



quadrupole degree of freedom as given in Ref. [pq| , where the spin operators J k are replaced 
by L fc = J fc + C k . 

In principle, the full kinetic term fl2"E| ) (after the reduction in terms of the parameter 
(l/o" 2 )) can be written down: it turns out still to be quadratic in the derivatives, with the 



coefficients given as power series in (/3/a). This is extensively discussed in Ref. |E5|. It 



is straightforward but cumbersome to obtain the corresponding coefficients of these higher 
terms in the kinetic energy. A much simpler procedure takes the higher contributions into 
account by redefining the coefficients of the lowest-order terms. In the example considered 
above, we would multiply the square bracket in Eq. (p8|) by a parameter (1/S 2 ), to be 
determined "experimentally" . This procedure is sensible as long as the expansion parameter 
(/? / cr) does not vary over a wide range . 

If we add to the terms in square brackets in Eq. ( [Z8D a contribution proportional to /3 2 , 
we obtain the Hamiltonian of a five-dimensional harmonic oscillator. In fact, such a term 
is available in (^), and we multiply it — following the above philosophy — by a factor C 2 
(there is another term relatively suppressed by a factor (1/cr 4 ), which we neglect). We can 
then write the whole expression as the number operator N2 of the quadrupole bosons times 
a constant depending on (C 2 /B 2 ), plus the zero-point energy. 

Taking into account all the terms in Eqs. (|28|) and (p9|) , the operator R 2 is given by a 
linear combination of the harmonic oscillator of the quadrupole degree of freedom discussed 
above, the eigenvalue C 2 of the Casimir operator of the SU(3)' group, the C 2 operator of the 
5*0(3)' group contained in S77(3)', a term proportional to a A and a constant. Therefore, as 
an ansatz for R 2 we propose 



R 2 = J^a 2 N 2 + k 1 (cr 4 - 36) + k 2 C 2 (\, fi) + k 3 ^2 C 2 k + h -] . (31) 
V ^2 k=1 



For convenience we have absorbed the zero-point energy of the harmonic oscillator and a 
term (36/ci + 1/4), in the parameter k± multiplying the constant term. 

Now we turn to the potential term V e ff in ([14]). In general, V e ff will be a complicated 

function of the three fundamental tensors l\ = g 2 , I 2 = ^1/5 j3 2 and ^3 = —^2/35 j3 3 cos(3c) 
(see Ref. [0), and also of those angles (p a b which are not excluded by gauge invariance. In 
keeping with the expansion in powers of (/?/cr) , we approximate V e jf by a function of g 2 
only. We also neglect any dependence on the angles (p a b which has to be considered as an 
ad hoc ansatz in the present context. There is, however, some motivation for this omission 
from the perturbative expressions calculated in Ref. [Q. We will further approximate the 
potential in the vicinity of its minimum by a linear combination of the terms 1/g 2 and g 2 , 
where the former can be assimilated into R 2 via a redefinition of the constant in Eq. fl3~T|). 
In fact, all the contributions to V e // which happen to have the same structure as one of the 
terms in R 2 , can be taken into account by redefining the corresponding parameters. The 
g 2 term in the potential can be interpreted as a gluon mass term as is clear from Eq. (|T9|). 
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We remark that, despite this mass term, the number of degrees of freedom of the spatially 
constant gluon field corresponds to transverse gluons due to the remnant of Gauss' law ([IT 
We then arrive at the final form of our ansatz for the effective Hamiltonian, 



H, 



1 d 2 



R 2 



eff 



2B Q Oq 2 



2B^ 



(32) 



with R 2 as in Eq. (|3T| ) and a new parameter C e from the modelling of the effective potential 
term. One of the main reasons behind the specific approximations made for the kinetic 
and effective potential terms is, of course, that the resulting Hamiltonian (^) is separable 
and analytically solvable. R 2 will be diagonalized simultaneously with H e ff, hence we can 
replace R 2 by its eigenvalue i? 2 , where \ refers to all quantum numbers which appear in 
this eigenvalue. H e ff then reduces to the Hamiltonian of a harmonic oscillator in g with 
continuous angular momentum, and the spectrum is 



k b = j|K 2 ^ + v + 4 



(33) 



with n e = 0, 1, 2, . . . and 



(34) 



Plugging in the eigenvalues for the different operators appearing in R 2 in Eq. ([IT]), we 
obtain the more explicit expression 



E 



n g ,N 2 




4\ —a 2 N 2 + &i(co>i + uj 2 + ou 3 )(ui + u 2 + uj 3 + 12) 

-D2 



+ k 2 {\ 2 + A/i + /i 2 + 3A + 3/i) + k 3 C(C + 1) + h 



4 f , 



(35) 



N 2 = and Uk = 
I). We see explicitly 



where N 2 is the quantum number of the five-dimensional harmonic oscillator, hence N 2 = 
0, 1,2,.. .. As we are interested only in energy differences, we have subtracted in ( |3~5"1) the 
eigenvalue for the vacuum, which has the quantum numbers n Q = 0, 
(k = 1,2,3), consequently A = /i = C = from Eqs. (fj), (|3|) and (| 
from (|35| ) that k 4 is constrained to take non-negative values. 

In our ansatz we have introduced a total of eight independent parameters which have to 
be adjusted in order that H e ff in (|32| ) approximate as closely as possible the exact effective 
Hamiltonian. Lacking knowledge of the latter, we will use information from lattice simu- 
lations for the same purpose. In the next section, we will determine the values of the six 
parameter combinations 

J-^-, h, k 2 , k 3 , hi, (36) 

V -Dp V -D2 

appearing in Eq. (p5|), by fitting the spectrum of H e ff to the glueball spectrum as measured 
in quenched lattice QCD. 
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Finally, we will give explicit expressions for the eigenstates of H e ff. Following Ref. ||12|| . 
we have 



\n e N2,AtL,JM,QC,(u 1 u 2 uj 3 )5) = 

^(e)^-A)/jWEWEW^ I ^0^U*)^$o)(*) . (37) 

K> K,K 

T^ e are the eigenfunctions of the one- dimensional harmonic oscillator in g given explicitly 
by 



KXq) 



2(B e C e ) 2+ »' 2 n 



{ng + v + 4) j ^ +?/2 < +3) ( v^> 2 ) eV^«V j (3 , 

where Lffi is a Laguerre polynomial and v was defined in ([34]). The function F in (|37j ) is 
the radial part of the five-dimensional harmonic oscillator and $ is the part that depends on 
c. Both are given in Ref. [p3]] . -D^ M represents the Wigner D-function [^(J whilst the other 
D-function is the representation matrix of 0(8) in the approximation of large a 2 . A has 
the meaning of seniority for the five-dimensional harmonic oscillator and gives the number 
of J = 2 gluon pairs not coupled to spin zero. The angular momentum L is restricted to 
the values allowed for a five-dimensional harmonic oscillator El, which in particular has no 



L = 1 state. The parameter t is a multiplicity index appearing in the classification of the 
five-dimensional oscillator and Q plays the role of the K in Eqs. (Q) and (|J) for the reduction 
of SU(3)' to S0(3)'. Finally, 5 is a further multiplicity index in the reduction of 0(8) to 
SU(3) colour. 

The basis given in Eq. (|37| ) is one of the possible realizations of the general expression 
given in Eq. @. In principle the direct relation could be found by using Dragt's theorem 
|2"8|| which relates the representation of states in terms of boson creation operators to the 
coordinate representation. We can in particular identify excitations of the one- dimensional 
oscillator (in g) and the five- dimensional oscillator (in (3, c) with gluon pairs created by the 
operators in (p. The relation is given through Eqs. ([T7|)-(|2TD and involves the parameters 
B g , C g and B 2 , C 2 , respectively. The total angular momentum of the pairs with quadrupole 
momentum is given by L (the other pairs do not contribute), where the allowed values of L 
correspond to the fact that these gluon pairs represent identical bosons. The minimum weight 
states in Eq. (||) transforming according to a certain 0(8) irrep (c^io^sO) , are described in 
([37]) by the corresponding 0(8) representation matrix. We can interprete the value of L as the 
angular momentum of these states, as is evident from considering a state with L = 0, where 
J = C by the Clebsch-Gordan coefficient in (|37D . In general, the total angular momentum 
J of a state arises as the (angular momentum) sum of L stemming from the J = 2 pairs and 
C originating from the minimum weight state. Interpreting the dependence of the energy in 
(P5|) on the quantum numbers of the states, the tensor gluon pairs and the unpaired gluons 
(described by the minimum weight states) interact strongly, while the scalar gluon pairs 
appear in our approximation as free particles. 

Using the basis (|37D, we can calculate the expectation value of g 2 which turns out to be 

<£ 2 > = ^-El , (39) 
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with E% taken from (j33l), i.e. including the vacuum energy. Eq. shows that the energy 
of a state is determined exclusively by the square of the intensity of the gluon field. Let us 
emphasize again that our model permits us only to calculate the spatial field distribution 
for zero-momentum states of glueballs, therefore we have no information about the internal 
structure of the states. 

We are confident that the Hamiltonian presented in ([32]) can serve as a first approximation 
to the exact effective Hamiltonian. In future developments the basis (|3"7|) will then play 
the important role of a perturbative basis in which to diagonalize the improved effective 
Hamiltonian. However, as there is at present no quantitative way to link our model to 
fundamental pure QCD it cannot be regarded as more than an ansatz somewhat in the spirit 
of a Landau- Ginzburg expansion of the free energy in statistical physics. We would claim 
that it is more concretely rooted in the fundamental microscopic theory than the former as 
our order parameter is directly related to the fundamental field content of QCD whereas, for 
example in superconductivity, the Landau-Ginzburg ansatz does not relate the scalar order 
parameter to the underlying electron degrees of freedom. 



4 The Glueball Spectrum 

The values of the parameters in Eq. (|35|) will be determined by adjusting the spectrum of 
our model to the masses of the 16 glueball states given by the lattice calculation of Ref. ||. 
For the three lowest-lying states we use the values cited in Ref. |29| (rather than Ref. ||) 
obtained as weighted averages over several lattice measurements done by different groups 



see the references in [2J|). The glueball state comes with a huge error in Ref. ||, so we 



decided to take the value from Ref. [BOfl instead, where a lattice Hamiltonian method was 



used. Finally, we use some of the (preliminary) values of Ref. |T0[ for the fit, which result 



from a new generation of glueball measurements on anisotropic lattices (see also Ref. ||31|| ' 



We have employed the CERN-MINUIT fitting routine [32| in order to obtain the best overall 



agreement. The resulting fit is given by the following formula (cf. Eq. (|35|)) 

E { nST ),C = °- 805 GeV x { 2n e + [i- 1 ^ 2 ^ + 0.165(0;! + u 2 + w 3 )(wi + u 2 + u 3 + 12) 
+ 0.877(A 2 + A/i + fi 2 + 3A + 3/x) + 0.028£(£ + 1)] 3 J . (40) 

Observe that the best value for the parameter is zero, so effectively we have only used 
five parameters. Note furthermore that the dependence on J arises only indirectly through 
N 2 , £ and the restriction given by the SU{2) Clebsch-Gordan coefficients in Eq. (j3?l). 

In Table 2 we list the results of the fit for the states considered in Ref. |J . For states 
marked with a cross only upper limits have been given in Ref. |J. Before commenting on the 
quantitative aspects of the fits, however, let us emphasize that the overall ordering of glueball 
states is to a large extent independent of the choice of the model Hamiltonian and relies only 
on the assumption that the number of constituent gluons is approximately a good quantum 
number and that more gluons correspond to higher energy. Of course, it is not clear a priori 
that the concept of a constituent gluon is a sensible one. However, in our approximation 
the quantity N = hi + h 2 + h 3 = 2{n Q + A^) + uJi + ^2 + ^3 is a good quantum number 
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and can be interpreted as the number of constituent gluons. It should be appreciated that 
the level ordering of the glueball spectrum does not follow any easily recognizable pattern, 
as was already stressed in Ref. ||. Using the above "minimal" assumptions, however, we 
already obtain good agreement with the ordering observed on the lattice. The two states 
that come out relatively too low, 1 and 3 , will be shifted upwards by large corrections 
stemming from the Casimir C 2 (\/-0 of the SU(3)' subgroup of the {7(3)' group related to 
0(8) (discussed in the previous section). Further below we will compare our results with 
those of other effective models, such as the MIT bag or the flux-tube model. 

Let us now proceed to a more quantitative comparison of our results with the lattice 
measurements. It must be noted that the overall agreement is surprisingly good, considering 
the simple-mindedness of our model. The most significant deviation of our values from the 
results of Ref. || occurs for the +_ state. However, the statistical error of the lattice 
measurement is probably too large for this discrepancy to be taken seriously. Furthermore, 
our value is in excellent agreement with the more recent calculation of Ref. ||10|| . As far as 



the latter is concerned, the statistical errors given in Ref. |I| are remarkably small, so if 
they are taken seriously, deviations from our values are more clearcut in this case. There is, 
however, an uncertainty in the overall energy scale for the measurements (not included in 
the errors cited here), so there might be a proportionality factor involved when going from 
one set of lattice measurements to another. Comparing our results to those of Ref. [K| we 
find appreciable deviations for the 3 +_ and 1 states, but most importantly for the 
where the deviation is twice as large as in all other cases. Our value is however in excellent 
agreement with the central value given in Ref. |3(J (after normalizing to the mass of the ++ 
according to Ref. P?|). 

Of course, deviations from our results are to be expected, since our model can only be 
considered as a first approximation to the exact effective Hamiltonian obtained from QCD 
by integrating out all non-constant modes. In particular, the form of the potential term 
(see the discussion preceding Eq. (P^D) is merely guessed at. As a consequence, one expects 
the degeneracy of several states in our model to be lifted by the correponding corrections. 
The good qualitative agreement with lattice results, however, leads one to believe that those 
corrections are not so large as to render our approximation meaningless. 

In Table 3 we have gathered the predictions from our model for states that may be 
accessible to future lattice measurements. Included are three excited states that have already 



been measured on the lattice fllQl , the results being in good agreement with our predictions. 
On a (hyper) cubic lattice one has to consider that rotational symmetry is broken down to 
invariance under the cubic group O. There exist only five irreps under spatial "rotations", 
termed A\, A2, E, T\ and T2. Consequently, the identification of a state on the lattice with 
a state in the continuum having angular momentum J is ambiguous and relies partially 
on the assumption that glueballs with higher J are heavier than those with smaller J, an 
assumption that is obviously not always fulfilled. Tables listing the O irreps contained in a 
given 5*0(3) irrep are available for example in Ref. |33| . 

State-of-the-art lattice calculations (on anisotropic lattices with improved actions) are 
able to determine the energies of the ground states and some first excited states for any 
given O irrep [[TTJ. Table 3 gives a list of our predictions for the energies of the first and 
second excited states for all R PC (R an irrep of O). Due to the degeneracy of several states 
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in our model, to be lifted by the inclusion of corrections to our ansatz for the effective 
Hamiltonian, the degree of excitation of a certain state (especially on the lattice) cannot 
always be determined from our mass estimates. In fact, these corrections together with mass 
shifts due to possible mixing of different states with the same 

jPC 

may change the sequence 

of states substantially in energy regions with a high density of states (see Figures 2-4). Our 
predictions for the level ordering in these regions should be taken with a grain of salt. 

We would like to point out that within our model we obtain a 4 ++ state lying considerably 
lower in mass than the 1 ++ . Given that both J = 1 and J = 4 contain an irrep T\ of the 
cubic group, on the lattice the 4 ++ would have been mistaken, had it indeed been found, 
for a 1 ++ according to the above-mentioned rule of thumb. Hence our value differs by more 
than two standard deviations from the result of Ref. ||, the latter being consistent with the 
more recent calculation in Ref. [|i~0[| . 

The states that are marked by an asterisk in Table 3 are only accessible on the lattice 
by measuring second excited states in the corresponding R PC sector, a technical challenge 
which has as yet not been met. More easily accessible states, to be compared with our pre- 
dictions, are T-f + *, A^~*, *, A%~*, E * and T{ * corresponding to the continuum states 
4 h , 3 H *, 2 + ~*, 3 *, 2 * and 1 *, respectively, with the predicted massses given in Table 
3. We should also mention that there are particular problems with the mass determination 
of ++ states on anisotropic lattices, which is why we did not cite the corresponding results 
in Tables 2 and 3. Ignoring these difficulties, the values given by the collaboration of Refs. 
0, [3TJ are (1.63 ± .03) GeV for the ++ and (2.84 ± .04) GeV for the ++ *. While the first 



value is consistent with the one given in Table 2 from Ref. |j29| , the second one is in good 
agreement with our prediction. 

In Fig. 2 we graphically compare our fit with the values for the glueball masses as cited in 
Refs. [^, [30]. The latter are represented by grey bars including the one-sigma deviation 
for the statistical errors. Light grey bars denote upper limits to the masses only. The other 
figures show our predictions for the full glueball spectrum up to a certain energy, obtained 
by using the parameters determined before. The states are ordered with respect to their PC 
eigenvalues and the 0(8) irreps. 

The general prediction then is that the spectrum is much denser than seen on the lattice, 
in particular in the sector PC = ++. Notice, however, that there is no contradiction with 
the lattice results due to the fact that on the lattice, with few exceptions, only the lowest- 
lying state in a given representation R PC is measured. A potential problem is the low-lying 
4 ++ as discussed above. 

Let us now discuss the present limits in energy for our prediction of the spectrum, repre- 
sented by the light grey bars in Figs. 2-4. While the parity of any state is simply given by 
P = (—1)^1+^2+^3 acC ording to Eq. (|7|), we are at present lacking a systematic determination 
of the charge conjugation eigenvalues C for states with more than six unpaired gluons, i.e. 
uji + uj 2 + 1^3 > 6 (there is work in progress on this matter ) . We will consider first the 
sectors with positive parity P = +1 and thus states with an even number of constituent 
gluons. Above 5.08 Gev, there are possible states with at least eight unpaired constituent 
gluons, the charge conjugation eigenvalues C of which are as yet undetermined. We do know, 
however, that the charge conjugation of states with 0(8) irreps (3320) and (4220) is C = — 1 
(see Appendix B). The lowest-lying of the O-undetermined states has (uiUJ2Uj 3 0) = (4330), 
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i.e. 10 unpaired constituent gluons, (A, fi) = (1,0), £ = 1 and n g = N 2 = 0, leading to 
J p = 1 + and a mass of 5.08 GeV. This state sets the limit indicated in Figs. 2 and 3. 

Turning now to the negative-parity sectors, there are states with at least seven unpaired 
constituent gluons, whose charge conjugation eigenvalues have yet to be determined. The 
lowest-lying among them has (uiUJ2UJ 3 0) = (3310), (A,yLt) = (0,2), C = and n e = N 2 = 0, 
and hence J p = 0~ and a mass of 4.46 GeV. All states with nine unpaired constituent 
gluons have masses above that value. We have, however, determined the charge conjugation 
of (3220) to be C = +1 and that of (4210) (not indicated in the figure) to be C = -1. Using 
the latter fact one can show that the lowest-lying 3 _+ state has a mass of 5.03 Gev, although 
there may be states with PC = — h and other angular momenta between 4.46 and 5.03 GeV 
that we cannot classify completely at present. However, none of these states can contain a 
representation T\ or A 2 on the lattice (i.e. J ^ 1,3,4, . . .). Similarly, we can conclude that 
the first state in this energy range containing A 2 ~ is the 3 ** at 4.68 GeV. However, there 
are candidate states for A±~*, A±~** and E ** below this value, so for the time being we 
cannot make predictions for these states on the lattice. All these facts have been used in 
constructing Tables 2 and 3. 

Finally, we would like to compare our results to those of other effective models |l|]-0. The 
oldest and also the most successful among these is the MIT bag model [TJ. All predictions 
for glueball masses in this model that we were able to find in the literature come out far too 
low. If we, however, take the predictions of Donoghue, Johnson, Li (and of Jaffe, Johnson 
for the (TM) 2 modes) [|IJ and simply scale the energy up by a factor of about 2.2, the lattice 
results are roughly reproduced. At this level, fine-structure corrections are absent, and the 
quality of the approximation is comparable to our model when we just order the states in 
accordance with the number of constituent gluons. 

Multiplying the energies by a factor of 2.2 roughly corresponds to changing B l / A to 
2.2 where B is the bag constant that characterizes the confining force. In the old bag 
model calculations this constant has been adjusted to the meson spectrum. It then appears 
plausible that B 1 ^ has to be changed by a factor of 2.2 « 9/4 in order to describe glueballs, 
since the factor 9/4 can be understood as the ratio of the colour charge of gluons to that 
of quarks. Of course, the rather large fine-structure corrections would have to be calculated 
anew to see if the MIT bag model can really account for a quantitative description of the 
glueball spectrum. As a first step, however, the masses of higher states should be estimated 
in order to make sure that the part of the spectrum above ~ 4 GeV as measured on the 
lattice can also be qualitatively understood within this model. The calculations in Ref. ]T5 



can be taken as a hint that one will have to consider states with more than three gluons. 

Let us remark that the constituent gluons of the bag model are very different from the 
ones in our model as they basically arise by putting free gluons into the bag. Consequently, 
boundary conditions at the bag surface and possible orbital motions play an important role. 
None of this makes any sense in our effective Hamiltonian description. A potential advantage 
of the bag model is that it allows, at least in principle, for a description of the inner structure 
of the glueballs. 

In Ref. |34j] it was suggested to consider the canonical dimension of the gauge-invariant 
operators that can create a glueball with a given J PC out of the vacuum in order to obtain 
a qualitative understanding of the level ordering of the spectrum. The results are somewhat 
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similar to the bag model (without fine structure corrections), but compare less favourably 
with the lattice data. Of course, it is also unclear how these ideas could be refined to yield 
quantitative predictions for the spectrum. 

There are other models like the "glue- lego" model @, the flux-tube model || and the 
non-relativistic potential model j|] which are, at least in their original forms, more or less 
excluded by the available lattice data. A relativistic version of the potential model is given 
by a recent Bethe-Salpeter calculation with massive constituent gluons [H. With a gluon 
mass of 0.4 GeV very reasonable values for the masses of the ++ ,2 ++ and h glueballs 
are obtained, while the 2 h comes out far too low and may be spurious. In Ref. |J a 
phenomenological QCD Hamiltonian is combined with a BCS type vacuum, which leads to 
an effective mass for the constituent gluons. Results are given for the masses of the ++ , 
2 ++ , _+ and 2~ + glueballs including excited states. The agreement with lattice data is 
satisfactory considering the approximations made. We emphasize that again the constituent 
gluons in all these models are different in nature from the ones in our model, and there is at 
present no way to relate the two concepts. 

Finally, sum rule calculations [B3] have their place somewhere between effective models 
and microscopic QCD. Up to now, predictions have been obtained only for the three lowest- 
lying glueballs and these are consistent with lattice measurements. 

In some of the phenomenological models a 1 _+ state with a relatively low mass is pre- 
dicted. This state is of some importance for experimental glueball searches, since it cannot 
be obtained from the quark model for mesons. In our model the 1~ + is a seven-gluon state 
and consequently its mass is rather large, in agreement with the lattice data. The problem 
with this state in other models can be traced to additional longitudinal degrees of freedom 
for the massive gluons, which in our case are absent as we have already pointed out in section 
3. 

We would like to conclude this section by remarking that in our model states with n e > 1 
(with the exception of the ++ glueball) appear as unbound states of n e ++ glueballs and 
the rest, insofar as the total energy of the state is just the sum of the energies of these 
components. Within the model, this is of course merely a consequence of the fact that 
the ^-dependence of the effective Hamiltonian is given by a harmonic oscillator. In physical 
terms, such states may correspond to very loosely bound states of n e ++ glueballs and other 
glueballs. On the lattice, these states may possibly be identified by analyzing their overlaps 
with different lattice operators and finite-volume effects. Such an enterprise is currently 



being envisaged [ID]. From our model, there is precisely one candidate for such a weakly 
bound state that has already been measured on the lattice, namely the h * state which 
should appear as a bound state of a ++ glueball with a h glueball. However, it must be 



noted that the lattice calculations of Ref. [10] give a larger mass for the + glueball and a 



slightly smaller one for the , so that the weak binding might be merely an artefact of 
our approximation. 

5 Conclusions 

In this paper we have performed a group-theoretic analysis of many-gluon states, based on 
the assumption that for the low-lying part of the spectrum only one mode of the gluonic 
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field (transforming under colour rotations and parity transformations like a constant mode) 
is important. The classification of the states involves a U(8) colour group and a U(3) spin 
group. The former is reduced via an 0(8) group to the colour group SU(3). For each 
given parity and charge conjugation the irreps of this 0(8) group classify bands of states. 
We remarked that in the context of charge conjugation the decomposition into elementary 
permissible diagrams (epds) proved very useful. The systematic determination of epds for 
higher 0(8) irreps is presently under investigation p^fl . 

One of the outcomes of the analysis is the possibility of defining a certain quantum num- 
ber which we identified with the number of constituent gluons, which again we emphasize are 
different from perturbative gluons or constituent gluons arising in other phenomenological 
models. The ordering of the states with respect to this quantum number allows for a quali- 
tative explanation of the level ordering of glueball states with respect to spin as observed in 
lattice calculations, assuming that the energy of the states increases monotonically with the 
number of constituent gluons. The ordering of the glueball states has been one of the most 
puzzling results of lattice simulations. It is worth noting that the MIT bag model might also 
be able to reproduce the glueball spectrum, provided that the value of the bag constant is 
chosen differently from the one used for quark states. To check this assertion new bag model 
calculations for glueballs are needed. 

In the second part of the paper we motivated a specific form for an effective Hamiltonian 
of low-energy pure QCD which depended only on the spatially constant modes of the gluon 
field. We introduced collective coordinates related to the monopole and quadrupole degrees 
of freedom of the intensity of the gluon field and rewrote the kinetic part of the Hamiltonian 



in terms of these variables using the results of Ref. [12]. Expanding in terms of b and 
(1/cr 2 ) where b parametrizes the absolute deformation of the intensity distribution and a 1 is 
related to the quantum numbers of the 0(8) group, and making an additional simplifying 
assumption for the potential term, we were able to deduce to lowest order a Hamiltonian 
which is a function of a one- and a five- dimensional harmonic oscillator and of Casimir 
operators contained in the U(8) D 0(8) D SU(3) chain of groups. Taking into account 
higher contributions via a redefinition of the coefficients of these operators we arrived at a 
phenomenological, QCD-motivated model Hamiltonian. The eigenstates of this Hamiltonian 
could be classified by the quantum numbers arising in the group-theoretical analysis, and in 
particular the number of constituent gluons turned out to be a good quantum number. 

By fitting the parameters appearing in the Hamiltonian we adjusted the spectrum of 
the model Hamiltonian to the "experimental" glueball spectrum, "experimental" in this 
context meaning lattice simulations. The general level of agreement was seen to be very 
good and predictions for several further states could be made. In qualitative terms, we 
predicted a much denser spectrum than that seen on the lattice without, however, running 
into inconsistencies with present lattice measurements. We have thus, with little input, 
principally the kinematical, i.e. group theoretical, structure of QCD, been able to deduce a 
great deal about the structure of an effective QCD Hamiltonian. 

We can naturally envisage several possible further investigations. Of course, an obvi- 
ous challenge is the derivation of an effective Hamiltonian from first principles. Without 
some form of infrared cutoff, such as a fermi-sized box, this will require a solution of the 
full crossover problem, i.e. a description of the interpolation between quarks and gluons as 
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effective degrees of freedom at high energies and hadrons at low energies. The collective 
coordinates introduced in this paper may prove useful in this context for the description of 
the low-energy degrees of freedom. Furthermore, a theoretic justification of the basic as- 
sumption made, that an effective Hamiltonian depending only on the constant modes of the 
field can yield the full glueball spectrum, is highly desirable. One could also think of possible 
extensions of the model, for example by including composite operators, in order to describe 
the inner structure of the physical states. Finally, it would of course be of interest to try and 
define collective variables related to quarks and anti-quarks, in a similar fashion to what has 
been done here for gluons. In that way the model could be extended to describe quark-model 
mesons and baryons, as well as the mixing between quarkic mesons and glueballs. Work on 
these topics is in progress. 
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Appendix 

A Reduction of 0(8) to SU(S) 

In this section will determine the reduction of the group 0(8) to SU (3). The reduction from 
£7(8) to SU(3) is known [|l|, as is the reduction from U(8) to 0(8) jl7j . The reduction 



from 0(8) to SU(3) can then be obtained in a recursive manner, starting from the simplest 
irrep [0] in U(8). 

For the scalar irrep [0] of U(8), the reduction to 0(8) is (0000) and the reduction to 
577(3) is (0,0), hence the S't/(3)-irrep (0,0) is contained once in the (0000) irrep of 0(8), 
which is trivial. Also the corresponding reduction of U(3) via SUj(3) to S*0(3) is trivial, i.e. 
the [0] irrep of U(3) only contains the spin J = 0. We have used the fact that the Young 
tableaux of the U (8) and U(3) irreps have to be the same. 

The U(8) irrep of one gluon (not presented in Table 1 because we give only those irreps of 
£7(8) which contain at least one scalar irrep of SU(3)) reduces to (1000) of 0(8) |l7j and to 
(1, 1) of SU(3), consequently the (1000) irrep of 0(8) must contain the (1, 1) irrep of SU(3). 
For the corresponding reduction of U(3) to S*0(3) we obtain J = 1, the expected result for 
a one-gluon state. 

For the case of a two-gluon system we have the symmetric [2] and antisymmetric [l 2 ] 
irreps of U(8). The first reduces to SU(3) as (2, 2) + (1, 1) + (0, 0) (see Ref. JKJ, |TJ), and for 
the corresponding reduction of U(3) to SO(3) we find J = 0,2, using the rules of Eq. (|3|). 
The reduction to 0(8) for [2] of 17(8) is given by (2000) + (0000) (Ref. ]TJ). The irrep (0000) 
already appeared in the reduction of the [0] irrep of U(8), and we know that it contains one 
colour zero irrep of SU(3). The other SU(3) irreps appearing in the list of the reduction of 
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U(8) to 577(3) are therefore contained in the (2000) irrep of 0(8), and in particular the latter 
cannot contain any colour scalar. For the antisymmetric irrep [l 2 ] of U (8) the reduction to 
SU(3) yields (1, 1) + (3,0) + (0,3) according to Ref. [|IJ, [L9], and the reduction to 0(8) is 
(1100), hence (1100) contains precisely the 577(3) irreps given above. 

In this manner, we can proceed recursively towards higher U (8) irreps, thereby obtaining 
the complete reduction of 0(8) to S77(3), and in particular the number of colour singlets 
contained in a given 0(8) irrep. The results for up to six constituent gluons are listed in 
Table 1, together with the possible values of J from the reduction to 50 (3). 



B Elementary Permissible Diagrams 

As claimed in section 2, any state characterized by the quantum numbers in (j^) can be 
obtained by a coupling of epds. Here we list all the epds with up to six gluons: 



(2, [2] 


,Ji) 


(3, [3] 


,Ji) 


(3, [I 3 ] 


,Js) 


(4, [2 2 ] 


,Ja) 


(4, [21 2 ] 


,Js) 


(5, [31 2 ] 


,Je) 


(5, [2 2 1] 


,Ji) 


(6, [3 2 ] 





dabc 
fabc 



Jl 



xbl 



bi 



[Ja] 
■h 
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[JA 



Y %1 

Y[3,l 2 } * faa'a"fabcd a 'b'c' 

Y[2 2 ,l\ * faa'a"fabcd a 'b'c' 
Y[ 3 2] * d, 



03 fai bi ci fa2 62 C2 fas 63 C3 

x bl F' 



[b\ x fet]W x [6 t x 6 tj 
Jx6t]W x[6 t,xfet,] 



x 61 



xbl 



■h 

w 
■h 



x [bl x bl 2 ]M x [bt 3 



[Js] 
Js 



(41) 



On the right-hand side the Young operator Y appears, which projects to a definite symmetry 
with respect to C7 (3) (or U(8)) and hence via [hih 2 h 3 ] = [o^o^s] to a definite 0(8) irrep. 
The spins denoted by (k = 1,2,..., 8) are limited to the values J\ = 0,2, J 2 = 1,3, 
J 3 = 0, J 4 = 0,2, J 5 = 1, J 6 = 0,2, J 7 = 1 and J 8 = 1,3. For the last three epds in Eq. 
(p|) intermediate couplings with spin labels J b and J c appear, for each of which one of the 
possible values has to be chosen. These are examples of the ambiguity in the choice of epds 
mentioned in section 2. 

The epds have been determined by a procedure described in Appendix D using the 
reduction of U(8) via 0(8) to 577(3) as deduced by the methods of Appendix A. For the 
explicit construction of the epds we have used the fact that every colour singlet state can 
be built up from the three "fundamental tensors" 5 a b, f a bc and d a t, c P5[ . As a result, the 
decomposition of all states with up to six consituent gluons in epds is known, and the 
charge conjugations of the states can be obtained easily from the charge conjugations of the 
corresponding epds. The latter are readily determined from the explicit expressions (HTf) as 
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shown in Appendix C. The results for the charge conjugations of all states with up to six 
gluons are given in Table 1, together with their parity eigenvalues determined by Eq. (|7]). 

In order to reproduce the lowest 1 _+ and H glueball states we have to consider minimum 
weight states with seven and eight gluons, respectively. The 1 _+ state is contained in the 
U(8) irrep [32 2 ], which reduces to the 0(8) irreps (3220), (3000), (2200), and others that 
do not contain colour singlets. Among the former, (3000) and (2200) can be obtained via 
a product of already determined epds coupled to a definite U(8) irrep. The (3220) irrep is 
given by a new epd, 



(7, [32 2 ],J=r 



^[S^ 2 ] * /aaia2^aa3a4^aibici^a262C2^ct3f)3C3 



4 x bi 



\ih 



[bl 



b 2 x b\ 2 



\[J2 



X 



'i >< b l 



xbl 



V34] 



[i] 



(42) 



which has parity P = — 1 and charge conjugation O = +1 (see Appendix C). 

The +_ state can be obtained by multiplying the two epds (3, [1 3 ],0) and (5, [31 2 ],0) 
coupled to the U(8) irrep [42 2 ] and the 0(8) irrep (4220). The charge conjugation of the 
state is the product of the charge conjugations of each epd, hence P = +1 and C = — 1. In 
view of the later discussion of excited glueball states, we note that we have determined in 
similar ways the parity and charge conjugation of colour zero states in the 0(8) irrep (4210), 

contained in the U(8) irrep [421], to be PC = , and for the 0(8) irrep (3320) contained 

in [3 2 2] we have obtained PC — H — . 



C Charge conjugation 

In order to determine the charge conjugation C of a state contained in a given 0(8) irrep 
(c^ic^^O) (with (0,0) as colour for its SU(3) subgroup) we introduced in section 2 the 
concept of epds. A list of epds for up to six gluons is given in (f4~T|) . Having determined 
the decomposition of a state into epds, its charge conjugation can be obtained simply by 
multiplying the charge conjugations of the epds. 

In order to deduce the properties of the epds under charge conjugation, we use the results 



in Ref. 21 . Consider the monomial 



T fe^V 2 * 6 int (A3) 

± a 1 a 2 fl„ w si u a 2 ■■■ u a n > ) 

where summation over the indices a*, is understood. Under charge conjugation the monomial 
transforms to 

V ai Va 2 -VaT aia2 aj>£b£...b£, (44) 

the values of the r\ a being given in Eq. @. 

As an example let us consider a special monomial of order four with colour zero and a 
certain angular momentum coupling, 

«,[Kx^x[6ixt!f{ ( . (45) 

Applying the charge conjugation operator and inserting a trivial factor tf a = 1, the SU(3) 
coupling dabcfade transforms to 

(VaVbVcdabc){VaVdVefade) = d' ab J' ade . (46) 
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Using Eq. we note that 

ji _ _j 

u abc ~ U'abc 

fade = fade > (47) 

i.e. the <i-symbol contributes a minus sign while the /-symbol contributes a positive sign. 

This can be readily extended to any monomial with colour zero. All one has to do is count 
the number of ci-symbols involved. Denoting this number by n^, the charge conjugation is 
given by 

C=(-l) nd . (48) 
The charge conjugations of the epds in ( p]) follow immediately. 



D Construction of States 

In this appendix, we will describe how states with up to six constituent gluons can be 
constructed from the epds listed in fl4ll) . At the same time, we will show how to actually 
determine the epds with the help of the reduction of Z7(8) via 0(8) to SU(3) as given in 
Table 1. 

We begin with the decomposition of a given state in a minimum weight state and gluon 
pairs. As stated in section 2, the lowest U(S) irrep containing a given 0(8) irrep (cj^^O) 
is described by the Young tableau [wiu^u^]. According to Eq. (||) we can build up all higher 
U (8) irreps with the same 0(8) irrep by applying an arbitary coupling of (ni + n 2 + 113) pair 
operators to the lowest state. Each pair is represented by the U(8) irrep [2], and is identical 
with the simplest epd (2, [2], J\) in the case M\ = J\ (maximum weight with respect to 
SO '(3)). 

A number of these pairs can be coupled to states represented by Young diagrams with 
an even number of boxes in each row. This restriction arises because the pairs themselves 
represent identical bosons with six degrees of freedom each. The problem of coupling of the 



pairs is equivalent to the reduction [17 



Z7(6) D 17(3) 
[1] - [2] 

[n] -> £[2m,2n2,2n 3 ], (49) 



where in the general reduction in the last line of Eq. ( BW ) the sum is over different partitions 
of n, n = n\ + n,2 + n^. For example, we have for the coupling of two and three pairs 

[2]x[2] = [4] + [2 2 ] 
[2] x [2] x [2] = [6] + [42] + [2 3 ] , (50) 

respectively. The resulting irreps then have to be coupled with a minimum weight state, 
[hih 2 h 3 ] = [uxu 2 uJ3\- 

We will now construct all the states appearing in Table 1 in this manner. Let us start by 
considering the (0000) irrep of 0(8). Coupling with one gluon pair operator, we obviously 
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obtain the U(8) irrep [2]. The results of the coupling of two and three pairs are given in Eq. 
([501). We thus get all the states with up to six gluons for the scalar irrep (0000) of 0(8) (cf. 
Table 1). 

In Table 1 we find two 17(8) irreps with three gluons, [3] and [l 3 ], both of which cannot 
be constructed by a coupling of pair operators. They hence give rise to epds, (3, [3], J 2) and 
(3, [l 3 ], J3), which — in contradistiction to (2, [2], Ji) — are minimum weight states. We 
now couple pair operators with these epds. As we exclusively consider states with up to six 
gluons, only one gluon pair will be coupled. We thus arrive at the U{8) irreps [5] + [41] + [32] 
for the 0(8) irrep (3000) and [31 2 ] for (1110) in agreement with Table 1. 

Similarly, there appear states with four and five gluons in Table 1, namely the U (8) irreps 
[2 2 ], [21 2 ], [2 2 1] and [31 2 ], which obviously cannot be obtained through coupling of the epds 
with two or three gluons, and hence give rise to new epds themselves (see Eq. (0)). We can 
still couple one gluon pair with the four-gluon epds to obtain the states given in Table 1. 

Finally, for the six-gluon states we consider couplings of the epds with three gluons as 
for example [3] x [3] = [6] + [42], where similar restrictions apply as for the coupling of 
several [2] irreps of 17(8) (see Eq. (106) of Ref. |]17|j ). We hence obtain new minimum weight 
states for the 0(8) irreps (6000) and (4200). The other couplings are [3] x [l 3 ] = [41 2 ] and 
[l 3 ] x [l 3 ] = [2 3 ], with the corresponding minimum weight states (4110) and (2220). 

There is one state in Table 1 which cannot be constructed in this way, namely the U(8) 
irrep [3 2 ], which consequently yields the last epd in (|4]). The decomposition of the states 
in epds is useful because it immediately determines the charge conjugation of the states in 
terms of the charge conjugation of the epds (see the previous Appendix C). 
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Table Captions 



Table 1: Reduction of U(S) to 0(8) up to six constituent gluons taking into account only 
those U(8) irreps which contain at least one colour scalar (the actual number of colour 
scalars is indicated in the third column). Also the reduction of U(3) via SUj(3) to 
the angular momentum group 5*0(3) is given, and we have indicated the correspond- 
ing values of J. Subindices refer to multiplicities which can be distinguished by the 
multiplicity index K. P and C refer to parity and charge conjugation, respectively. 

Table 2: Glueball masses calculated within the model compared to lattice calculations. 
The values of the parameters used in our model are given in Eq. (f5|). See sections 2 
and 3 for our classification of the glueball states displayed in the first four columns. 
Columns five and six give the usual J PC classification along with the value of the mass 
determined from our model Hamiltonian. 

The lattice results which have been used to fit the parameters are given in the last 
two columns. While the next to last column cites the best confirmed values currently 
available, the last column refers to preliminary results of a new generation of glueball 
calculations on anisotropic lattices (we did not include the value for the ++ glueball 
due to problems with this specific state in the new method). The errors displayed refer 
exclusively to statistical errors, while there might be an additional overall variation 
in the scale. In particular, the scales used in the last two columns may slightly differ 
from one another. Values which are only upper limits, are marked by a cross. Some of 
the lattice values have actually been read off the figures in the respective publications, 
hence the data given might be erroneous in the last digit. 

Table 3: Predictions from our model for states that can possibly be measured on the lattice 
in the near future. The first five columns are as in Table 2, while the sixth gives 
the corresponding representations on the lattice, where the angular momentum J is 
replaced by the label R referring to the representations of the cubic group O (see Ref. 
[j33|l). A superscript asterisk marks an excited state in columns five and six (in column 
six we have listed only those components which can be seen on the lattice as at most 
second excitations). Where the degree of excitation of a certain state could not be 
determined from our mass estimates due to the degeneracy of several states in our 
model, possible higher degrees of excitation were indicated by superscript asterisks in 
parenthesis. 

The last column gives our mass estimates and a few preliminary results from the 
lattice. An asterisk in front of a value indicates that the corresponding state can only 
be determined from a second excitation on the lattice, which appears to be rather 
demanding with present techniques. 
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Figure Captions 



Figure 1: Comparison of our predictions and current "best values" for the glueball masses 
as given in Table 2 (columns six and seven), arranged according to their PC eigenvalues. 
Our results are indicated by solid lines, while the lattice values, including statistical 
errors (one-sigma deviations), are represented by grey bars. Light grey bars indicate 
upper limits. More than two-sigma deviations occur for the mass of the 1 ++ versus 
our prediction of the 4 ++ (see the discussion in section 4) and for the H , where our 
value is however in excellent agreement with the data of Ref. [[10] (see Table 2). 

Figure 2: Our prediction for the complete glueball spectrum with PC = ++ below 5.08 
GeV, using the parameters in fl4*0|). Note how dense the spectrum is in comparison 
to what one would see on a lattice (Table 3). States listed in Table 2 and Figure 1 
are marked by bold face letters. The lowest-lying state is, of course, the vacuum (not 
indicated in the figure) which has {ujiUj^oj^) = (0000), (A,//) = (0,0), C = and 
n e = N 2 = and hence J PC = ++ and whose energy has been set to zero. 

Above 5.08 Gev, there are possible states with at least 8 unpaired constituent gluons 
{uj\ + co> 2 + co>3 > 8). The charge conjugation eigenvalues C of these states are as yet 
unknown. 

Figure 3: As in Figure 2, but for PC = H — . Note that the states with (wiu^u^O) = (3320) 
and (4220) have 8 unpaired constituent gluons. In these particular cases we have 
determined the charge conjugation to be C — — 1. 

Figure 4: As in Figure 2, for states with PC = — h and and masses below 4.46 GeV 

(with the exception of the (cjiu^sO) = (1110) states at 5.03 GeV, see section 4). 
Between these two values, there are states with at least 7 unpaired constituent gluons, 
whose charge conjugation eigenvalues have yet to be determined. We do know, however, 
the charge conjugation of states with (uiUiurfS) = (3220) to be C = +1. 
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